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ABSTRACT 

We present ultra-high resolution spectra for a set of nearby F-G-K stars on, or close to, the main 
sequence. The wavelength shifts of stellar lines relative to their laboratory wavelengths are measured 
for more than a thousand Fe i lines per star, finding a clear correlation with line depth. The observed 
patterns are interpreted as convective blue-shifts that become more prominent for weaker lines, which 
are formed in deeper atmospheric layers. A morphological sequence with spectral type or effective 
temperature is apparent. Two K giant stars have also been studied. The velocity span between weak 
and strong lines for these stars is larger than for the dwarfs and subgiants of similar spectral types. 
Our results show that convective wavelength shifts may seriously compromise the accuracy of absolute 
spectroscopic radial velocities, but that an empirical correction may be applied to measured velocities. 



L INTRODUCTION 

Convection in late-type stars penetrates into the pho- 
tosphere, producing inhomogeneities. Temperature varia- 
tions of up to a few hundred Kelvin are apparent in optical 
images of the Sun, and a velocity field of several kilometers 
per second is observed in spatially resolved solar spectra. 
The warmer upflows appear as bright granules surrounded 
by narrower, cooler downflows. While present technology 
cannot resolve stellar granulation, there is evidence of its 
presence in all late-type stars. The temperature and veloc- 
ity inhomogeneities produce absorption line profiles which 
are shifted and asymmetric, effects that are readily noticed 
in ultra-high dispersion stellar spectra. 

An extensive literature exists on solar line asymmetries 
and shifts (see, e.g., Neckel & Labs 1990; Asplund et al. 
2000; Pierce & Lopresto 2000). Modern stellar studies 
were pioneered by Dravins (1974, 1985) and Gray (1980, 
1981). Dravins (1999) provides a short review of recent 
work. Line asymmetries and line shifts are consequences 
of the same phenomenon and, therefore, give largely equiv- 
alent information. When the number of lines measured is 
limited, but the available features are mostly unblended, 
line asymmetries carry more information. On the contrary, 
if the spectrum is heavily crowded, line shifts, which can 
be measured for many lines, are likely to provide more 
information. A limited wavelength coverage, largely dic- 
tated by the detector size and the scarcity of high-accuracy 
laboratory wavelengths, steered stellar studies toward the 
analysis of line asymmetries. 

Line bisectors - a measure of the line asymmetry - vary 
smoothly across the HR diagram (Gray & Toner 1986; 
Dravins 1987; Gray & Nagel 1989; Dravins & Nordlund 
1990). Other parameters, such as rotation, chemical com- 
position, magnetic fields, or binarity, are likely to play a 
role, but the limited data available have prevented an in- 
depth study. A recent comprehensive work on the Fc i 
spectrum (Nave et al. 1994) has provided accurate lab- 
oratory wavelengths for 9501 lines. These data and im- 
provements in astronomical instrumentation make it fea- 



sible and practical to turn to line shifts as a complement 
to line bisectors. 

Gravitational and convective wavelength shifts system- 
atically affect absolute determinations of radial velocities. 
Other effects can also introduce systematic radial veloc- 
ity errors, but they are expected to be less important for 
most stars (Lindegren, Dravins, & Madsen 1999). Grav- 
itational redshifts are proportional to the mass-to-radius 
ratio. For spectral lines formed in the photospheres of 
dwarf stars, the gravitational shift is the same for all lines, 
in the range between 0.4 and 2 km s~^. When accurate 
parallaxes and photometry are available, comparison with 
evolutionary models allows us to determine dwarf masses 
and radii within 8 % and 6 % , respectively (Allende Pri- 
eto & Lambert 1999). It is therefore possible to estimate 
the gravitational shift within 10 %, yielding worst-case 
uncertainties of ^ 0.2 km s~^. Convective line shifts may 
vary with spectral type. The velocity difference between 
weak and strong lines is about ~ 0.6 km s~^ for the Sun 
(Allende Prieto & Garcia Lopez 1998) and probably more 
than 1 km s~^ for the F5 subgiant a CMi (Procyon; Al- 
lende Prieto et al. 2002), but no systematic study across 
the HR diagram has been published. Obviously, accurate 
spectroscopic studies of absolute radial velocities cannot 
afford to neglect these shifts. 

In this paper, we analyze optical spectra of several late- 
type stars obtained with the High Resolution Spectrograph 
(HRS) coupled to the Hobby-Eberly Telescope (HET). 
The high quality and large spectral coverage of the spectra 
allow us to measure line shifts for a large number of Fe i 
lines. Section |^ describes the observations and the data. 
In Section ^ we describe the analysis, and in §^ we discuss 
the results. 

2. OBSERVATIONS 

Our observations were obtained with the 9.2m HET at 
McDonald Observatory and the HRS (TuU 1998). The 
HRS saw first light on March 17, 2001, starting science 
observations shortly afterwards. This fiber-coupled echelle 
spectrograph uses an R-4 echelle, cross-dispersing gratings. 
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we determine the horizontal scatter between the polyno- 
mial fit F(A) and the observed spectrum /(A). For each 
observed wavelength Ai considered in the fit, we find nu- 
merically the roots of the polynomial F{X) — f{Xi) using 
Laguerre's method, select the appropriate root, and then 
derive the rms scatter between the observed and fit wave- 
lengths. Assuming the errors are Gaussian, this scatter is 
divided by the square root of the number of points enter- 
ing the polynomial fit to derive the uncertainty as- As an 
example, application to the solar flux spectrum provides a 
nearly normal distribution of uncertainties that peaks at 
0.00025 A and shows bumps on the long-value tail, signal- 
ing blends. Nave et al. (1994) classify their wavelengths 
into four categories, depending on their uncertainty (crj), 
which ranges from 0.4 mA to more than 10 mA. Most lines 
have errors below 5 niA and about half of the lines have 
errors below 1 mA^. 

The uncertainties in the measured wavelength shifts, 
a — \/(t'^ + cr'f , are then converted to velocity, and their 
distribution n{a) is determined. Then, we use n{a) to pre- 
dict the distribution of errors (S) in the observed velocity 
shifts as 
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and a mosaic of two 2kx4k CCDs. In the selected mode, 
HRS provides a FWHM resolving power R ~ 120, 000 and 
complete spectral coverage between 4094 and 7890 A, with 
the exception of the gap (5977-5999 A) between the CCDs. 
The observations were performed in queue mode in May 
and June 2001 (see Table 1). An observation consisted 
of two exposures that were extracted and combined. The 
final signal-to- noise (S/N) per pixel was > 400 at most 
wavelengths. 

The spectra were processed with the tasks in the IRAF 
echelle package. We removed the bias signal, flatfielded 
and extracted the spectra, subtracting the low background 
in the wide interorder regions. In addition to the program 
stars observed with the HET, we have included in this 
study the spectra of the Sun (Kurucz et al. 1984) and a 
CMi (AUende Prieto et al. 2002). 

3. ANALYSIS 

3.1. Stellar parameters 

The effective temperature (Toff) has been estimated 
for all the program stars by one or more authors using 
the Infrared Flux Method (Blackwell & Lynas-Gray 1994; 
Alonso et al. 1996, 1999). The T^s of a Boo (Arcturus) 
has been independently derived by di Benedetto (1998) 
and Griffin & Lynas-Gray (1999). We have adopted the 
values in Table 1. All our stars were analyzed by AUende 
Prieto & Lambert (1999), and we accepted the surface 
gravities derived there from the (B — V) colors and the 
Hipparcos trigonometric parallaxes. The iron abundances 
in Table 1 are a straight average of previous spectroscopic 
analyses compiled by Cayrel de Strobel et al. (2001). A 
recent discussion of the preferred stellar parameters of a 
CMi can be found in AUende Prieto et al. (2002). Fei 
lines at 4602.00 and 4602.94 A with accurate <? /-values and 
damping parameters were analyzed with synthetic spectra 
computed from flux-constant LTE model atmospheres to 
estimate the projected rotational velocities in Table 1. 

3.2. Measuring the line centers 

A variety of procedures have been used previously to 
determine the central wavelength of a spectral line. The 
most popular methods involve some type of least-squares 
polynomial fitting or polynomial interpolation (see, e.g., 
Neckel & Labs 1990; AUende Prieto & Garcia Lopez 1998). 
As the stellar spectral lines are asymmetric, the details of 
the measuring procedure affect the derived central wave- 
length. Polynomial least-squares fits are well-suited for 
error estimation, however, there is some arbitrariness in 
the selection of the ideal order and wavelength interval to 
use in the fits. For the slowly rotating stars analyzed here, 
we have found that a reasonable choice, the same for the 
whole sample, provides stable results regardless of the ex- 
act selection for the order and wavelength interval. Our 
particular choice is to fit 35 mA around the line center 
with a third order polynomial. Preceding the fitting by a 
cubic-spline interpolation with a step of 5 mA proved use- 
ful for the HRS spectra, but degraded slightly the quality 
of the results for the superior solar atlas. 

To estimate the error in a derived central wavelength, 

^This comparison indicates that an extension and improvement of the laboratory wavelengths for neutral iron and other species is desirable. 
A spectrograph used to obtain stellar spectra may be well suited for this purpose (see AUende Prieto 2001) 
^With the exception of a CMi, for which only 588 lines were measured 



where C is chosen such that N{S)dS — 1. The de- 
rived N{S) can be compared to the observations in order 
to constrain the intrinsic scatter. In all cases, the derived 
and measured histograms are similar. 

3.3. Searching for correlations 

The number of iron lines identified and measured for 
each star varies from 1340 to 1569^. Solar analyses have 
shown clear correlations of the line shifts with line strength 
(AUende Prieto & Garcia Lopez 1998) and wavelength 
(Hamilton & Lester 1999). We explore here the variation 
of the shifts with line strength, as this is the dominant ef- 
fect. The line equivalent width and the central line depth 
have been the most widely used line strength indicators. 
As pointed out by Pierce & Lopresto (2000), the use of 
the equivalent width avoids saturation for the strongest 
lines. On the other hand, equivalent width measurements 
are distorted by blends in the line wings which may hardly 
affect line depths. Blends are a particularly serious prob- 
lem when dealing with cool stars, as some of those in our 
sample. Therefore, we have adopted the line depth. 

The left-hand panels in Fig. |l] show the line velocity 
shifts for the stars in our sample as a function of the resid- 
ual flux at the center of the line /c {— 1— line depth). In 
these panels, the shifts of the individual lines are marked 
with dots. Grouping the lines in bins of 0.1 in depth, and 
after applying a median filter, we derive the average val- 
ues shown as asterisks with error bars. The solid lines are 
3'^-order polynomial fits to guide the eye, and the bro- 
ken lines trace the original median values. We have not 
attempted to disentangle the stellar radial velocities from 
the gravitational shifts and the atmospheric motions. The 
absolute values on the vertical axis have been shifted to 
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Fig. 1. — The left-hand panels show the relative velocity shifts of Fe I lines with respect to their rest wavelengths. The velocity scale for 
each star is shifted to have a null velocity for the median shift of the strongest lines considered. The dots show the shifts for the individual 
lines. The asterisks with error bars are mean values after applying a median filter, and the solid line is a polynomial fit to them. The dashed 
curve shows the median shifts for the different bins. The right-hand panels show histograms of the scatter of the velocity shifts around the 
polynomial fit. The expected scatter (see text) is shown by the symmetric solid curve. The mean velocity shifts determined from the analysis 
of the a Boo atlas of Hinkle et al. (2000) are displayed as filled circles in the left-hand panel for this star. 



Table 1 
Observed sample 



Star Obs. date Sp. Type Tes'^ log 51 [Fc/H] vsmi S{v) 

(K) (dex) (dex) (km s"!) (km s"!) 



a CMi 


McDonald atlas 


F5 IV 


6530 


3.96 (0.02) 


-0.05 (0.03) 


2.8 (0.1) 


0.88 


Sun 


FTS atlas 


G2 V 


5770 


4.44 (0.01) 


0.00 


1.9 (0.1) 


0.54 


70 Vir 


June 5 


G4 V 


5455 


4.04 (0.11) 


-0.11'' 


3.2 (0.1) 


0.72 


/Lt Her A 


May 28 


G5 IV 


5523 


3.87 (0.04) 


0.23 (0.07) 


4.1 (0.2) 


0.59 


6 Aql 


June 11 


G8 IV 


5040 


3.56 (0.07) 


-0.02 (0.10) 


3.1 (0,3) 


0.57 


'/ C(>p 


May 29 


KO IV 


4944 


3.40 (0.18) 


-0.12 (0.14) 


3.1 (0.1) 


0.51 


a Boo 


June 9 


K1.5 III 


4290 


2.00 (0.56) 


-0.51 (0.03) 


3.0 (0.5) 


0.93 


C Dra 


May 29 


K2 III 


4430 


2.31 (0.32) 


-0.13 (0.04) 


5.1 (0.5) 


1.04 



"^Uncertainties in Tefj are 100 K or smaller 
single measurement is available 
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have a null median velocity shift for the strongest lines; 
a condition that the polynomial fits were as well forced 
to satisfy. The right-hand panels show the scatter of the 
line shifts about the average points. The scatter expected 
from Eq. (1), mainly the result of the uncertainties in the 
laboratory wavelengths, the finite width of the lines, and 
the presence of photometric noise, is also displayed (thick 
symmetric curve). The c^S) given in the right-hand pan- 
els corresponds to a Gaussian fit by least-squares to the 
observed scatter. The velocity span of the polynomial fits, 
S{v), is listed in Table 1. 

A spectroscopic atlas of a Boo at a superior resolving 
power and S/N was published by Hinkle et al. (2000). 
The corresponding left-hand panel in Figure |l| compares 
the average shifts from our measurements from the atlas 
(filled circles) with those from the HRS spectrum (aster- 
isks with error bars). The agreement is quite satisfactory. 
The a of the scatter in our spectrum and in the atlas agree 
within 10 m s~^. 

4. DISCUSSION 

The pattern of the velocity shifts is similar for all the 
stars, and in qualitative agreement with solar results. The 
net convective blue shifts of the lines strengthen toward 
deeper photospheric layers, and therefore affect the weak 
lines more than the strong ones. Despite our limited sam- 
ple, there is an indication of a smooth dependence of the 
average velocity shifts with spectral type. This effect is 
clearer in our analysis of line shifts than in previous stud- 
ies using bisectors. As argued in the introduction, this 
is likely the result of line shifts measurements being less 
affected by blending than line bisectors in late-G and K 
type stars. 

Our preliminary results indicate that the trend of line 
shifts as a function of line strength can be determined as 
a function of spectral type and gravity. At least for late- 
type dwarfs, assuming the strongest lines in the spectrum 
are free from convective shifts (as is the case for the so- 
lar photosphere), it is possible to correct for convective 
wavelength shifts to within ~ 0.2 km s^^. 

The measured range of wavelength shifts for neutral iron 
lines is larger in a CMi than in cooler stars of classes IV 
and V. The range of wavelength shifts is also larger, reach- 
ing up to ~ 1 km s~^, for the K giants than for the G-K 



dwarfs and subgiants in our sample, which is not surpris- 
ing. Although their lower effective temperatures provide 
less flux to transport, their lower atmospheric pressure im- 
plies larger convective velocities. Fig. 1 suggests that, 
once we account for the observational errors, there is little 
room for intrinsic line-to-line scatter in the F-G-K stars 
sampled here. Convection is supposed to cease in main 
sequence stars with spectral types earlier than about F2, 
but other velocity fields must be present in those atmo- 
spheres, judging from the asymmetry of the spectral lines 
(e.g. Gray & Nagel 1989). Determining accurate absolute 
radial velocities demands an understanding of the wave- 
length shifts in the spectra of these stars too. 

Solar observations have shown differences in the center- 
to-limb variation of the granulation along the central 
meridian and the equator (e.g. Beckers & Taylor 1980; 
Rodriguez Hidalgo, CoUados, & Vazquez 1992). Observa- 
tions of a large number of stars may then show peculiarities 
at a given spectral type depending on the orientation of 
the rotational axis. Stars with enhanced magnetic fields 
are also expected to be peculiar in terms of observed wave- 
length shifts, as the enhanced magnetic fields may hinder 
the convective motions. Furthermore, line shifts may vary 
with time for stars exhibiting an analog of the solar 11 
year cycle. 

Observations of line shifts for a significant number of 
stars with high quality are required for a deeper under- 
standing of granulation in stellar atmospheres, its relation- 
ship with other stellar phenomena, and the role of surface 
convection in the structure and evolution of stars. Our first 
results show that measurements of convective line shifts 
are a must in order to derive accurate absolute radial ve- 
locities. 

We thank the HET staff for their outstanding job mak- 
ing possible science observations with HRS since day one. 
The Hobby-Eberly Telescope is operated by McDonald 
Observatory on behalf of The University of Texas at 
Austin, the Pennsylvania State University, Stanford Uni- 
versity, Ludwig-Maximilians-Universitat Miinchen, and 
Georg-August-Universitat Gottingen. NSO/Kitt Peak 
FTS data used here were produced by NSF/NOAO. This 
research was supported in part by the NSF (grant AST- 
0086321). 
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